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Foreword 


The  effects  of  the  oceanic  environment  on  acoustic  propagation, 
sensor  systems,  and  weaponry  have  become  widely  recognized  as 
important  in  prosecuting  modem  naval  warfare.  Cost  and  time  are 
prohibitive  for  conducting  field  experiments  to  measure  these 
effects  under  all  the  conditions  that  should  be  expected.  Simulation 
models  of  the  oceanic  environment  can  produce  realistic  data 
approximations  that  can  be  used  in  designing  and  testing  systems  and 
tactics.  This  report  describes  the  first  ocean  simulation  model 
developed  at  the  Naval  Oceanographic  and  Atmospheric  Research 
Laboratory  and  provides  potential  users  with  the  details  necessary 
to  implement  and  operate  the  model. 


W.  B.  Moseley  (j.  B.  Tu^z,  Captain,  USN 

Technical  Director  Commanding  Officer 


Executive  Summary 
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Oceanographic  simulation  models  are  valuable  tools  for  imitating  tempera¬ 
ture,  salinity,  density,  sound  velocity,  and  other  environmental  areas  of  naval 
interest.  The  Naval  Oceanographic  and  Atmospheric  Research  Laboratory 
was  tasked  to  implement  and  operate  ocean  simulation  models.  This  manu^ 
presents  the  theoretical  foundations  for  the  first  version  of  the  model,  the 
details  of  the  primary  model's  structure,  and  the  computer  programs  and 
instructions  for  implementing  and  using  the  model. 

The  first  version  of  the  model  was  based  on  the  premise  that  the  most 
important  features  of  interest  to  the  Navy  are  produced  by  internal  waves. 
Subsequent  versions  of  the  model  will  incorporate  variability  produced  by 
shear,  fronts,  mixed-layer  processes  and  turbulent  mixing.  A  fuU  three- 
dimensional,  time-varying  model  is  expected  in  about  3  years. 
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Ocean  Simulation  Model  for  Internal  Waves 


I.  Introduction 

A.  Background 

Over  the  past  two  decades,  acousticians  have  real- 
I  ized  that  acoustic  propagation  is  greatly  affected  by 
the  variability  in  the  internal  wave  field  and  by  the 
variability  of  currents  and  sound  velocity  on  the  larger 
scales  of  fronts  and  eddies,  (Rattd  and  Tappert,  1975; 
Dvson  et  al.,  1976;  Dcsaubies,  1976a,  1976b,  1978; 
liwart,  1980,  1989;  Macaskill  and  Ewart,  1984).  This 
I  significant  discovery  led  to  the  need  to  obtain  simu¬ 
lations  of  oceanic  density,  current,  and  sound  velocity 
fields.  One  method  for  studying  the  acoustic 
fluctuations  relies  on  a  direct  simulation  of  the  ocean 
environment.  A  statistical  model  of  the  environment 
can  also  be  directly  incorporated  into  an  acoustic 

•  model. 

Simulations  are  used  to  model  the  performance  of 
proposed  instruments,  often  before  these  instruments 
are  designed  or  constructed.  This  type  of  modeling  is 
particularly  useful  when  the  instrument  senses  or 
^  processes  the  signal  nonlinearly.  Use  of  such  simu¬ 

lations  can  save  time  and  money  in  the  design  and 
construction  of  sensor  systems. 

Oceanographic  simulations  also  have  great  potential 
utility  in  designing  sampling  strategies  for  hydro- 
graphic  surveys.  A  field  of  simulated  data  can  be 
P  easily  constructed  and  one  sampling  method  tested 
against  another.  Proper  application  of  this  type  of 
testing,  if  conducted  prior  to  a  cruise,  could  save 
several  days  of  ship  time  per  cruise. 

An  ocean  simulation  model  can  be  defined  as  a 
method  (an  algorithm  implemented  on  a  computer) 

*  that  will  provide  a  representation  of  a  set  of  fields  of 
oceanic  variables  over  a  given  space-time  domain  at 
a  specified  resolution.  The  algorithm  can  simply  be  a 
method  for  searching  a  data  base,  or  it  can  employ 
sophisticated  methods  for  generating  stochastic  fields 

^  with  specified  statistics. 

'  B.  Project  Objectives 

The  Naval  Oceanographic  and  Atmospheric 
Research  Laboratory  (NOARL)  was  tasked  to  develop 


and  test  a  sequence  of  simulation  models  for  the  ocean 
simulation  project.  The  end  product  will  be  the  capa- 
bihty  to  simulate  temperature,  salinity,  density,  sound 
velocity,  and  current  fields  in  any  ocean  region  over 
any  reasonable  time  and  space  domain  and  resolution. 
The  simulations  are  intended  for  use  by  the  acoustics 
and  oceanographic  community  to  provide  a  well- 
defined  testbed  for  evaluating  environmental  effects 
on  acoustic  systems,  and  by  the  weapons  design, 
nonacoustic  ASW  and  oceanographic  communities  for 
similar  purposes. 

The  first  version  of  the  simulation  model  is  based 
on  the  premise  that  the  features  important  to  the 
Navy  are  produced  by  internal  waves.  These  waves 
typically  have  periods  between  15  minutes  and 
24  hours  and  horizontal  wavelengths  on  the  order  of 
tens  of  meters  to  tens  of  kilometers.  The  waves 
vertieally  advect  constant  density  surfaces  to  produce 
time-  and  space-varying  fields  of  temperature, 
salinity,  and  sound  speed. 

Subsequent  versions  of  the  simulation  model  will 
incorporate  variability  produced  by  shear,  fronts, 
mixed  layer  processes,  and  turbulent  mixing.  A  full 
three-dimensional,  time-varying  model  is  expected  in 
about  3  years. 

C.  Purpose  of  this  Manual 

This  manual  provides  the  simulation  model  user 
with  information  that  will  assist  in  understanding  the 
theoretical  and  numerical  bases  for  the  simulation, 
the  details  of  program  structure  needed  for  imple¬ 
mentation  and  maintenance,  and  provides  instructions 
for  using  the  program  suite  to  generate  simulations. 
This  manual  describes  the  algorithms  used  in  the 
simulation  and  discusses  some  of  the  less  obvious 
numerical  methods.  The  details  of  the  program  struc¬ 
ture  are  presented  and  a  guide  to  program  operation 
follows.  The  program  codes  are  listed  in  NOARL 
Technical  Note  No.  059.  ASCII  source  code  can  be 
obtained  from  NOARL ’s  Physical  Oceanography 
Branch  on  5V4-inch  diskette  or  V2-inch  magnetic  tape. 
For  copies,  please  write  to:  Commanding  Officer, 
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Attn:  Code  331,  Naval  Oceanographic  and 
Atmospheric  Research  Laboratory,  Stennis  Space 
Center,  MS  39529-5004. 

II.  Background  and  Derivation 
of  Algorithms 

A.  Stochastic  Representation  of  Internal 
Wave  Field 

Internal  Waves.  The  internal  wave  field  is  com¬ 
posed  of  motions  in  the  frequency  band  between  the 
inertial  frequency,  f,  and  the  Brunt-Vaisala  frequency 
N.  (Strictl>  speaking,  this  is  the  band  in  which  freely 
propagating  internal  waves  are  permitted — forced  or 
evanescent  internal  waves  can  exist  outside  this 
frequency  band.)  The  sobriquet,  internal  waves,  also 
implies  the  existence  of  a  dispersion  relation,  relating 
wave  frequency  to  wavenumber. 

The  best  known  (and  most  widely  used)  model  of 
the  internal  wave  field  was  first  formulated  by  Garrett 
and  Munk  (1972).  This  model  is  based  on  an 
exponential  approximation  to  the  Brunt-Vaisala 
frequency  profile  and  an  approximate  fit  to  observed 
spectra  and  cross  spectra  of  density  and  velocity  fields. 

Tides.  After  internal  waves,  the  tides  are  the  second 
most  important  sources  of  rapid  fluctuations  in  the 
ocean.  They  can  be  decomposed  into  a  barotropic 
and  baroclinic  component.  These  have  widely 
separated  length  scales.  The  barotropic  tide  has  an 
extremely  large  length  scale  when  compared  to  the 
internal  (baroclinic)  tide. 

Barotropic  Tide.  The  barotropic  tide  has  such  long 
horizontal  length  scales,  uiat  for  most  of  our  modeling 
effort,  it  may  be  taken  as  constant  or  slowly  varying 
horizontally  and  uniform  vertically.  The  tidal  current 
vectors  trace  out  an  ellipse  in  the  horizontal  plane 
over  a  tidal  cycle.  The  primary  problem  in  modeling 
the  barotropic  tide  is  to  choose  the  amplitude  of  the 
motion.  For  the  open  ocean,  the  amplitude  is  taken  to 
be  about  1  m,  and  currents  of  about  10  cm/sec  are 
assumed. 

The  barotropic  tide  is  also  important,  as  it  interacts 
with  the  bottom  to  generate  an  internal  tide  (of  the 
same  frequency  as  the  barotropic  tide)  and  higher 
frequency  internal  waves,  such  as  lee  waves  observed 
at  Stellwagen  Bank  in  Massachusetts  Bay  and  propa¬ 
gating  solitons  in  the  Western  Mediterranean  Sea. 

Internal  Tide.  The  internal  tide  is  thought  to  be 
generated  by  the  coupling  of  the  internal  wave  field 
with  the  barotropic  tide  in  regions  where  the  bottom 
has  a  very  sharp  discontinuity,  such  as  the  continental 
shelf  break.  The  internal  tides  formed  in  these  regions 


are  radiated  into  the  deep  ocean  as  beams  of  internal 
wave  energy  at  tidal  frequencies.  They  often  proceed 
shoreward  as  internal  solitons  under  favorable  density 
stratification.  The  generation  regions  are  sometimes 
found  at  midocean  banks  or  seamounts. 

The  largest  amplitudes  of  the  internal  tide  occur 
on  the  continental  shelves,  with  the  amplitudes  in  the 
open  ocean  being  much  smaller.  DeWitt  et  al.  (1986) 
have  studied  the  internal  M2  tide  in  the  Rockall 
Channel  area,  west  of  the  British  Isles.  They  found 
that  internal  tides  are  consistent  with  the  model 
of  Prinsenberg  and  Rattray  (1975),  with  most  of 
the  energy  being  in  the  first  vertical  mode.  The 
Prinsenberg  and  Rattray  model  indicates  that 
the  amplitude  of  the  nth  mode  should  be  distributed 
as  (1/rt)  ■  sin(n  71 //]///),  where  //j  is  the  depth  of  the 
shallow  region  and  H  is  the  deep-ocean  depth.  Local 
information  is  not  available  at  present,  so  the 
amplitudes  of  the  internal  tidal  modes  are  assumed  to 
be  proportional  to  the  inverse  mode  number. 

Outline  of  Basic  Garrett-Munk  Model.  In  this 
simulation  effort,  the  latest  version  of  the  “Garrett- 
Munk  spectrum”  is  used.  There  are  a  number  of 
different  versions  (Garrett  and  Munk,  1972,  1975, 
Cairns  and  Williams,  1975;  Desaubies,  1976b, 
Munk,  1978).  The  Munk  (1978)  formulation  is  used 
because  it  incorporates  more  observational  data  than 
the  others,  and  it  is  defined  in  a  frequency-mode 
formulation,  which  is  more  appropriate  to  the 
eigenmode  expansion  approach.  A  table  of  terms  and 
definitions  is  provided  for  convenience  (Table  1). 

The  following  assumptions  are  either  made 
explicitly  or  implicitly  in  the  derivation  of  the  Garrett- 
Munk  (GM)  model. 

•  The  Brunt-VSisala  frequency  profile  is  given  by 
Nq  exp(-z/t),  where  Nq  is  the  Brunt-Vaisala 
frequency  at  the  surface,  z  is  the  depth,  and  6  is  a  scale 
depth.  For  the  usual  GM  formulation,  Nq  =  3  eph  and 
b=  1300  m. 

•  The  fields  of  interest  are  given  by  a  sum  of 
weighted  wave  eigenfunctions  based  on  a  separated, 
linear,  internal  wave  equation  without  shear. 

•  Local  WKBJ-like'  scaling  is  valid  for  all  modes. 

•  The  spectra  are  consistent  with  observed  internal 
wave  spectra. 

•  The  spectrum  is  separable  by  mode  and  frequency. 

•  The  energy  associated  with  each  mode  can  be 
“smeared”  in  the  (o,  k  plane. 

The  spectral  representations  used  in  the  most  recent 
formulations  (Munk,  1978)  for  the  spectra  of  vertical 


'See  Niyfeh  (1973)  for  «  discussion  of  the  WKBJ  approximiuon  method. 


Table  1.  Terms  and  definitions. 


b 

Scale  depth  in  the  GM  model.  Usually  1.3  km. 

T 

horizontal  range. 

B(oj) 

GM  frequency  weighting  function. 

s„ 

Spectrum  of  horizontal  velocity. 

ev' 

Brunt-vaisaia. 

Spectrum  of  vertical  velocity. 

E 

Nondimensional  GM  energy  level,  6.3x10'*. 

s; 

Spectrum  of  vertical  displacement. 

1 

Inertial  frequency  -  2  =  sin  (latitude). 

^(Z) 

Weighting  function  used  in  the  Numerov  method. 

1(x) 

Function  of  x. 

U 

Horizontal,  eastward  velocity  component. 

Fu 

GM  form  of  the  horizontal  kinetic  energy  spectrum. 

u, 

Horizontal  u-velocity,  y'th  mode. 

F. 

GM  form  of  the  vertical  kinetic  energy  spectrum. 

V 

Horizontal,  northward  velocity  component 

GM 

Garrett-Munk. 

Horizontal,  v-velocity,  yth  mode. 

H 

Bottom  depth. 

GM  mode  weighting  function. 

w 

Vertical  velocity  component. 

HU) 

w, 

X 

Vertical  velocity,  yth  mode. 

Horizontal,  eastward  coordinate. 

HKE 

Horizontal  kinetic  energy 

y 

Horizontal,  northward  coordinate 

i 

Mode  number. 

2 

Vertical  coordinate. 

j* 

Mode  scale,  usually  3. 

a 

Direction  of  a  section. 

k 

magnitude  of  horizontal  wavenumber. 

Weight  of  yth  mode  at  frequency  co. 

m 

Vertical  wavenumber. 

Vertical  displacement. 

n 

Direction  index. 

K 

horizontal  wavenumber  vector. 

N(Z) 

Brunt-Vaisaia  frequency  at  depth  z. 

0„y(<o)  Random  direction  of  rrth  wave,  y  fh  mode,  frequency  co. 

No 

GM  base  Brunt- vaisaia  frequency,  usually  3  cph. 

CD 

Frequency. 

NH 

Vertically  integrated  Brunt-Vaisaia  frequency. 

(|»„y(a>)  Random  phase  function. 

velocity,  F^,  and  horizontal  velocity,  as  functions 
of  frequency,  O),  mode,  j,  and  local  Brunt-VaisSia 
frequency,  are 

Fw  ( CO,  j)  =  N-\z)  ( 0)2  -/2)  E  ( (0.  j)  ,  (1) 

(CO,  j)  =  b^N^Niz)  (0)2  +  /2)  o)-2£  (0),  j )  ,  (2) 


where 


E  ((0,y)  =  £B(ci))//0), 

(3) 

B(a))  =  27i  ‘/w''  (0)2  ^ 

(4) 

//(;■)=(/ 

(5) 

6=  1300  m, 

(6) 

No  =  3  cph  , 

(7) 

£  =  6.3x  10^ 

(8) 

The  mode  scale  j,  is  assumed  equal  to  3. 

Extensicn  to  Realistic  Brunt-Vaisala  Profiles.  To 
extend  the  GM  model  to  general  profiles  of  Brunt- 
Vaisaia  frequency,  a  number  of  further  assumptiorts 
have  been  made. 


•  The  displacement  and  velocity  fields  can  be 
expressed  as  a  sum  of  weighted  wavefunctions, 
integrated  over  the  following  frequencies;  from  f  to 
min^ge  (max,  [A/(z)]). 

•  The  wavefunctions  are  determined  by  the  sepa¬ 
rated  z-dependent,  linear,  internal  wave  equation  and 
local  Brunt- Vaisaia  frequency  profile,  N{z). 

•  The  vertically  integrated  energy  in  each  mode 
has  the  GM  form  in  frequency  and  mode. 

•  In  the  first-level  model,  there  are  no  shear  effects. 

These  assumptions  result  in  an  energy  scaling 

somewhat  different  from  the  GM  model.  This 
scaling  is  discussed  in  the  following  section. 

Energy  Scaling.  The  general  form  of  the  eigen¬ 
function  expansion  in  terms  of  the  wave  functions 
may  be  summarized  as: 

[N 

w  (z,  r,  a.  t)  =£2;  j  (to)  H^  (q);  z)  • 
j  «  Jf 

exp  {/  [xkjcos  Q„jia))  +  ykjSin 

+  (ot  +  ,  (9) 

where  x  =  r.  cos(a),  y  =  r .  sin(a),  W,  ((o;z)  =  vertical 
mode  of  order  j  at  frequency  to,  and  is  the  solution  to 
the  ordinary  differential  equation, 
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(10) 


integrated  over  the  entire  water  column,  can  be  written 
in  terms  of  the  scaled  horizontal  velocities  U  and  V: 


Wj  +  k^B[N(z),(O,f]Wj=0, 
Wj  =  0  at  z  =  0,H. 


The  primes  in  equation  10  denote  differentiation  with 
respect  to  z  and 


B[N{z).w,f] 


N^(z)-(ji'^ 


(II) 


where 

py(tD)  is  contribution  from  mth  mode  at  frequency  to; 

r  is  range; 

a  is  direction  of  section; 

0„^(to)  and  ^„j((o)  are  independent  random 
variables,  uniformly  distributed  over  [0,2n); 

Q„ji(o)  is  the  direction  in  which  the  nth  wave  of 
mode  j  is  propagating; 

(to)  are  random  phases  associated  with  specific 
modes,  directions  and  frequencies. 

Here,  we  employ  the  formalism  of  Pierson  (1955) 
to  represent  the  ensemble  of  realizations  of  a  random, 
stationary,  Gaussian  process.  (See  Pierson,  1955,  or 
Kinsman,  1965,  for  a  detailed  explanation  of  the 
Vdto  formalism.) 

The  primary  problem  in  reconciling  the  GM 
spectral  model  of  the  internal  wave  field  arises 
from  the  departure  of  real  ocean  density  profiles  from 
the  exponential  profile  assumed  by  Garrett  and  Munk. 
Because  real  density  profiles  have  scales  of  variation 
much  smaller  than  the  (approximate)  vertical  wave¬ 
length  of  a  particular  mode,  the  WKBJ-like 
approximations  used  in  the  GM  (Garrett  and 
Munk,  1972)  model  cannot  be  fully  justified.  Much 
of  the  formalism  of  the  GM  model  may  be  kept, 
however,  by  appropriate  scaling  of  the  wave  functions 
in  a  manner  consistent  with  the  vertically  integrated 
GM  spectra.  The  derivation  of  this  scaling  follows. 

Define  the  orthonormal  modes  of  the  vertical 
velocity  as  W-  and  subject  to  the  relation; 


f' 


W,Wj  [{N  (z)  -  0)  )/(0)  -/2)]  dz  =  ,  (12) 


HKE  = 


=  V2  (  Uj  +  Vj  dz 

Jo 


(13) 


1  i-  -^W.^dz  (14) 

0  0)2  ^2  dz  ' 


HKE  =  '/2 


r../i 

Jo 


w/  dz  (15) 


HKE  =  1/2  1  +  £  . 


to 


(16) 


From  the  GM  formulation,  equation  2  (Munk,  1978), 


f2 

nKE  =  b^N„  \+J—  £(to.;)  N(,z)  dz  .  (17) 

Jo 


0) 


However,  if 


NH=l  N{z)dz, 

h 


(18) 


then 

Py  (©)  =  Ib^N^NH  Eicaj)  (19) 

is  obtained  as  the  proper  scaling  for  the  vertical 
eigenfunctions  of  vertical  velocity. 

Here  b  =  1300  m,  Nq  =  3  cph  =  0.005236  rad/sec, 
H  =  30(X)  m  (for  the  test  case  presented  in  the  section 
on  usage)  and  E (co,y )  =  B((o)  H{j)  as  defined  in 
equations  3-5. 

If  all  the  random  phase  terms  are  combined  into 
V^fto),  then  the  expression  for  w  is 


where  W,  and  Wj  are  solutions  to  the  internal  wave 
equation  10  for  Afferent  modes. 

Then  define  the  scaled  {Uj,  Vj,  Wj  )  as  {Uj,  yj,Wj) 
=  Vj,  Wj).  The  specific  horizontal  kinetic 

energy  (HKE)  of  the  yth  mode  at  frequency  to. 


W  (z,  f)  =  £  w^(z,  t) 
] 


p^(co)VV/z;(ol 


exp  {i[(or-t  v/co)]}  Vdm  . 


(20) 
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The  autocorrelation  function  (Kinsman,  1965)  can 
be  constructed  as 

z)  z)> 


=  V2  Pj(co)Wj(2; 

Jo 


w)  cos[a)(r- /')]  ifo) ,  (21) 


where  (co)  is  set  to  zero,  where  (o  <  /  and  (£>>  N. 
The  spectrum  of  the  vertical  velocity  fluctuations  of 
the  yth  mode  at  a  given  level,  z,  are  thus  given  by 

z,j)=  pj (w)  Wj((a, z) ,  (22) 

and  the  total  vertical  velocity  spectrum  by  the  sum 
over  aU  modes  is 

5>,z)=lP)  H'2((o,2)  .  (23) 


The  GM  relations  to  the  spectra  of  displacement  and 
horizontal  velocity  are 


(0^  -  P  (P 


spectrum  is  computed  for  the  jth  mode.  From  inte¬ 
grating  Equation  1, 


VKE  =  2b^N^  {(P-P)  E(0) 


J)f 

Jo 


(z)  dz  (27) 


is  obtained.  But,  by  integrating  the  squared  vertical 
mode  velocities,  we  obtain 


VKE  =  >/2p 


r" 

7  (CO)  Wj 

Jo 


=  b^N^NHEim 


These  can  be  consistent  only  if 


,;■)[  wjdz. 

Jo 


i(p  -p ) 


Jo 


{,z)dz=NH  Wj{z,(i))  dz  .{29) 


From  the  orthogonality  condition. 


-f 


{(P-P)  =  {N^-(P)  Wjiz,  to)  dz  ,  (30) 


((o,z)  =  to'^  (to.z)  .  (25) 

If  ((o,z)  is  known,  then  integrating  it  in  the  vertical 
with  the  inner  product  weighting  function  gives 


lo  to'-/' 


5,^  (to,  z )  dz 


I  P'((o)  f 
^  Jo 


(0^-p  ' 


=  I  p'((0)  =  2fe2;V^/V//I  £(to,y) 

j  j 

=  2b^N^NHBmiHj 


=  2b^N^N  . 


we  have 


J(A^-(o2)VF2(z,to)dz  {  N-U 
0  Jo 

=  /V//  lFj(z,to)dz. 

Jo 


,  to)  dz  . 


However,  if  N{z)  is  set  to  equal  a  constant,  then  the 
equality  does  not  hold,  so  the  scaling  (26)  is  not  con¬ 
sistent  with  the  GM  spectrum  for  the  vertical  velocity. 

The  source  of  this  disagreement  can  be  traced  to 
the  WKBJ-like  approximation  for  the  vertical 
wavenumber,  m  =  kN(oP  - /  and  the  relation  that 


_  cP-p 


which  arises  from  the  relation  between  the  derivatives 
of  w  and  the  horizontal  components  of  velocity. 


There  is  a  logical  inconsistency  here,  however. 
Suppose  the  vertical  kinetic  energy  from  the  GM 
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TTius,  if  the  GM  form  for  the  vertically  integrated 
spectrum  of  the  HKE  is  taken  as  our  basis  for  scaling 
the  vertical  eigenfunctions,  then  a  spectrum  of  vertical 
kinetic  energy  inconsistent  with  the  GM  form  is 
obtained.  The  consistent  form  for  the  spectrum  of  the 
vertical  kinetic  energy  is  given  by  equation  23,  and 
the  spectrum  for  the  vertical  displacements  is  still 
related  to  by  equation  25.  The  spectrum  for  the 
HKE  can  be  obtained  from 

(0),  zj)  =  (1  +/2/C02)  wy  2 

and 

(0),  z ,)  =  X  (to,  z ,7 ).  (34) 

J 

B.  Plotting/Interpolation  Methods 

The  density  of  the  output  data  is  not  sufficiently 
high  to  produce  smooth  plots,  so  the  data  must  be 
interpolated  onto  a  finer  mesh,  consistent  with  the 
screen  resolution.  The  interpolation  is  performed  by 
a  bilinear  function  defined  over  each  rectangle  in  the 
original  output  mesh.  This  method  has  one  drawback: 
it  does  not  ensure  continuity  of  first  derivatives.  This 
flaw  is  compensated  for  by  the  increased  speed  and 
does  not  significantly  affect  the  visual  aspect  of  the 
plot  if  the  original  resolution  is  moderately  high. 


C.  Export  Algorithm 

For  plotting  the  fields  on  an  IBM-compatible  per¬ 
sonal  computer  (PC)  using  a  1 6-color  video  graphics 
adapter,  the  data  have  been  encoded  as  two-character 
numerical  ASCII  groups,  varying  from  00  to  99.  These 
groups  are  transmitted  over  the  local  area  network 
via  Kermit  or  other  convenient  file  transfer  protocol 
in  blocks  of  39  pairs  and  plotted  on  a  PC.  A  header  in 
the  file  allows  for  reading-in  scaling  and  plot 
generation  information.  A  typical  plotting  program 
for  a  PC  is  listed  in  NOARL  Technical  Note  59. 

A  more  efficient  coding  algorithm  is  possible,  as 
only  16  colors  are  required.  The  data  could  be  encoded 
for  color  level  on  the  generating  computer  and 
transmitted  with  two  color  codes  per  byte,  an 
improvement  of  a  factor  of  4  in  storage.  More  efficient 
packing  methods  are  possible  and  may  be  used  in  the 
future,  but  most  of  these  methods  require  much  more 
computadonal  time  at  both  the  source  and  receiver 
and  are  prone  to  transmission  errors. 


in.  Numerical  Methods 

A.  Eigenfunction  Calculations 

The  eigenvalues  and  eigenfunctions  for  the  internal 
wave  differential  equation  are  computed  using  the 
Numerov  method.  This  recursive  finite-difference 
method  has  the  advantage  of  being  very  fast  and  is 
fifth-order  accurate.  (The  routine  was  originally 
written  by  Dozier  (undated  manuscript)  of  Scientific 
Applications  International  Corporation  (SAIC)  and  is 
a  part  of  the  original  code  supplied  by  Rubenstein, 
also  of  SAIC.) 

The  Numerov  method  for  the  differential  equation 
y"+f(z)y  =  0  is  given  by  the  relation 

[1  T(z  +  Az)]y(z  -f-  Az)  +  [\  +  T{z-  Az)]y(z  -  Az) 
=  [2  10  T(z)]  y(z)  -  AzV‘^V240  -i-  .... ,  (35) 

where  r(z)  =  Az^  /(z)/12. 

The  integration  proceeds  from  the  lop  and  the  bot¬ 
tom  toward  a  point  zq  in  the  oscillatory  region  of  the 
solution  (i.e.,  where  /(z)  >  0).  The  eigenvalue  is 
updated  until  a  matched  solution  is  found  that  satisfies 
continuity  of  the  function  and  its  first  derivative  at 
the  matching  point  to  the  desired  degree  of  accuracy. 
The  mode  is  checked  by  ensuring  that  the  number  of 
zero  crossings  of  the  solution  matches  the  mode 
number.  This  method  occasionally  fails  when  the 
eigenvalues  are  extremely  close  together,  but  this 
problem  does  not  appear  to  be  serious  in  most 
circumstances. 

B.  Integration  Method 

Several  methods  are  available  to  estimate  integrals 
using  numerical  techniques.  The  question  arose  as  to 
which  would  provide  a  better  estimat  of  the  simu¬ 
lated  fields.  The  result  of  a  number  of  tests  was,  that 
due  to  the  random  phase  function,  no  convergence 
exists  other  than  a  convergence  to  the  Gaussian 
probability  distribution  function  over  an  ensemble  of 
realizations  of  the  integral.  We  found  that  the  simplest 
form  of  numerical  integration,  rectangular 
approximation,  gave  results  as  good  as  any  other 
method  and  was  more  efficient. 

C.  Choice  of  Integration  Subintervals 

If  a  uniform  grid  is  chosen  over  which  to  perform 
the  numerical  integration,  then  the  grid  must  be  made 
very  fine  to  adequately  resolve  the  structure  of  the 
GM  spectrum  near  the  inertial  frequency.  To  use  a 
coarser  mesh  while  still  resolving  the  spectrum, 
a  logarithmic  transformation  is  performed  on  the 
(0  axis.  This  ojxration  requires  that  the  functional  form 


computed  over  the  sum  be  modified  and  consideration 
taken  in  the  choice  of  the  differentia]  weights. 

The  transformation  is  given  by 

b  fb  fin  b 

fix)  dx^\  X  if)  d{ln  X)  =  ey  f(ey)  dy  . 

I  J a  J In  a  (36) 

The  form  on  the  left  is  exactly  equal  to  that  on  the 
right,  but  if  there  is  a  singularity  in  f(x)  at  x  =  0,  and 
a  is  near  0,  then  the  form  on  the  left  provides  a  greater 
accuracy  for  the  same  number  of  terms  in  the  sum  or, 
conversely,  less  terms  for  the  same  accuracy.  As  an 
example  for  fix)  =  x'^,  a  =  0.001,  b  =  5.0,  Gaussian 
quadrature  did  not  reproduce  the  answer  to  better  than 
30%,  while  the  form  on  the  right  agreed  with  the  true 
answer  to  within  the  roundoff  error  of  the  computer. 
This  transformation  is  easily  modified  for  singular¬ 
ities  in  f(x)  at  points  different  from  x  =  0.  For  the 
estimation  of  the  internal  wave  field,  the  logarithmic 
form  has  been  chosen,  because  the  singularities  in 
the  GM  spectra  are  located  at  (o=/«  A. 

When  computing  the  integral,  it  was  necessary  to 
compute  the  locations  of  the  ends  of  each  subinterval 
in  the  transformed  space  to  properly  resolve  the  square 
root  of  the  differentials. 

D.  Problem  of  Probability  Distribution  Due 
to  Finite  Summation 

The  problem  encountered  here  arises  from  the 
necessity  of  approximating  the  integral  (9)  as  a  finite 
sum  over  a  set  of  intervals  on  the  (o  axis.  The 
probability  distributions  of  each  component  of  the  sum 
are  Cauchy  distributions;  thus,  the  sum  will  be  a 
sum  of  Cauchy  distributions.  As  the  sum  approaches 
the  limit  of  the  integral,  the  probability  distribution 
of  the  sum  along  the  time  axis  or  across  the  ensemble 
at  given  times  will  approach  a  Gaussian  distribution 
almost  surely  (i.e.,  as  the  limit  in  the  mean). 

The  computation  has  been  modified  somewhat  to 
produce  a  rapid  approach  to  a  Gaussian  distribution 
with  few  terms  in  the  sum.  The  method  is  to  add  a 
small  random  phase  that  varies  linearly  in  time  to  the 
constant  (in  time)  random  phase-that  is,  to  replace 
0„^((o)  in  equation  9  with  <t>„/(o)  +  e  y„y(to)  t,  where 
e  is  a  small  number  and  and  y„y((o)  are  random 
variables,  uniformly  distributed  over  (0,  2n]. 

IV.  System  and  Program  Structure 

A.  System  Structure 

Simulation  Model.  The  design  philosophy  for  the 
simulation  model  is  based  on  simulating  processes  in 


the  ocean  (or,  alternatively,  regions  in  to,  k  space)  in 
order  of  their  energetic  importance  on  the  time  scales 
of  interest,  typically  scales  of  hours  to  days.  Because 
the  most  energetic  part  of  the  oceanic  spectrum 
appears  in  the  internal  wave  band,  modeling 
the  internal  wave  field  has  begun,  and  is  followed  by 
internal  tides,  mesoscale  features,  fronts  and  eddies, 
and  large-scale  features,  in  that  order.  The  larger  scale 
features  arc  important  tc  the  smaller  scales  in  that  they 
control  the  dynamics  of  the  internal  wave  field  through 
the  mean  density  field  and  mean  shear  field,  but  they 
will  be  incorporated  in  future  versions  of  the  model. 

The  first  version  of  the  simulation  model  was 
designed  to  provide  global  geographic  coverage  of 
the  ocean  and  seasonal  variations.  Several  databases 
exist,  but  the  Levitus  5°  database  (Levitus,  198'*''  was 
chosen  for  the  initial  modeling  based  on  storage 
considerations.  This  database  is  used  to  define  strati¬ 
fication  by  providing  spatially  dependent  profiles  of 
temperature  and  salinity  from  which  the  density  field 
is  computed.  The  oceanic  variability  is  produced 
by  advecting  the  temperature  and  salinity  fields  by 
vertical  wave  motions  at  given  positions  and  times, 
based  on  a  variation  of  the  GM  internal  wave  model. 

To  make  the  first-level  ocean  simulation  model 
generally  useful,  it  was  necessary  to  ensure  that  the 
model  could  be  successfully  run  on  a  variety  of 
computers  and  that  the  output  data  would  be 
compatible  between  different  machines  and  operating 
systems.  This  requirement  has  forced  restrictions  on 
the  source  code  to  ensure  portability  of  the  program 
and  compatibility  of  the  data.  At  present,  the  model 
has  been  run  on  VAXA'MS  and  CONVEXAJNIX 
systems.  The  restrictions  are  not  onerous,  however. 
They  consist  primarily  of  adhering  to  the  Fortran  77 
conventions.  Some  differences  still  exist  between 
operating  systems,  that  prevent  the  code  from  being 
completely  portable  (e.g.,  the  CONVEX  uses  the 
convention  that  the  record  length  parameter  in 
opening  direct  access  files  is  given  in  bytes,  whereas 
the  VAX  convention  specifies  the  record  length 
parameter  in  4-byte  words). 

Program  MODEL  1  is  version  1  of  the  simulation 
model.  Its  output  is  a  time  series  of  vertical  sections 
of  vertical  displacement,  temperature,  salinity,  and 
three  velocity  components.  A  diagnostic  file  provides 
timing  information  on  the  VAX  but  not  on  the 
CONVEX. 

Export  Programs.  Programs  are  available  for 
exporting  the  ocean  simulation  model  output  from 
binary  format  into  ASCII  format.  While  binary 
formatted  files  require  less  storage  space,  problems 
might  arise  when  transferring  the  files  to  computers 
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with  different  hardware  architectures.  For  example, 
when  transferring  files  from  a  VAXA'MS  system  to 
a  CONVEX  system,  the  byte  order  must  be  reversed. 
Binary  characters  can  also  be  interpreted  as  control 
characters.  In  ASCII  format,  the  data  are  more  easily 
transferred  to  other  computers  due  to  the  non¬ 
machine-specific  structure  of  ASCII  files. 

Program  MODEL  l_EXPORT  computes  model 
field  data  produced  by  the  ocean  simulation  model  for 
plotting  on  PCs.  Interactive  prompts  query  for  neces¬ 
sary  input  and  the  program  produces  an  ASCII  file. 
Plots  can  then  be  produced  on  the  PC  using  the  data 
in  the  ASCII  file. 

Program  MODEL  l_SV_EXPORT  uses  the  ocean 
simulation  model  binary  output  file  to  produce  an 
ASCII  file  containing  sound  velocity  data  for  transfer 
to  other  computers. 

Program  ASCLEV  converts  the  Levitus  data  file 
from  binary,  direct  access  format  into  ASCII, 
sequential  format. 

Program  LEVASC  converts  the  Levitus  data  file 
from  ASCII,  sequential  format  into  binary,  direct 
access  format.  Binary  format  is  used  by  the  oceaii 
simulation  model. 

Diagnostic  Programs.  Several  diagnostic  programs 
are  available  for  evaluating  model  output.  This  suite 
of  programs  provides  a  variety  of  functions,  which 
includes  generating  data  listings,  creating  plots,  and 
examining  energy  levels. 

Program  MODEL  l_PLOT  computes  and  plots 
model  fields  on  the  DEC  GKS  Workstation  2000.  No 
files  are  produced.  This  program  is  intended  to  provide 
a  high-resolution,  quick  look  at  the  data  but  suffers 
from  the  lack  of  a  hard-copy  device  on  the  GKS 
workstation. 

Program  MODEL  l_CONTOUR  computes  and 
plots  model  fields  as  contours.  A  DISSPLA  metafile 
is  produced  as  output.  The  metafile  may  be  plotted 
using  the  DISSPLA  postprocessor. 

Program  MODELl .ENERGY  computes  internal 
wave  energies  as  plot  fields.  Horizontal  and  vertical 
distributions  are  plotted.  DISSPLA  metafiles  are 
produced  as  output.  The  metafile  may  be  plotted  using 
the  DISSPLA  postprocessor,  DISSPOP. 

Program  MODEL  1_UVW_READ  reads  the  .AUX 
and  .UV  output  files  from  the  ocean  model.  The 
velocity  fields  U,  V,  and  W  are  displayed  on  the  user’s 
terminal  screen. 

Program  MODEL  1_CTL_LIST  lists  the  control 
information  file  MODELl.CTL.  This  file  con¬ 
tains  information  relating  to  program  restart  and 
eigenmode/eigenvalue  control  and  use. 


Program  MODEL  1  .LOOK  lists  model  displace¬ 
ment  fields  on  the  user’s  terminal  screen.  This 
program  is  primarily  intended  as  a  means  of 
dmermining  whether  the  typical  values  of  the  fields 
are  reasonable. 

Program  MODEL  l.EIGLOOK  lists  and  plots 
eigenmodes  and  eigenvalues  of  the  internal  wave- 
functions  produced  and  stored  by  MODEL  1.  The  plots 
are  produced  as  a  DISSPLA  metafile  and  can  be 
plotted  using  DISSPOP. 

Program  MODESUB.TEST  computes  the  internal 
wave  modes  for  a  Brunt-Vaisaia  frequency  profile 
and  plots  them.  The  plots  are  generated  as  DISSPLA 
metafiles  that  can  be  plotted  using  DISSSPOP. 

Program  DISPERSION  produces  modal  dispersion 
diagrams  for  a  Brunt-VSisSla  frequency  profile  at  a 
given  location  and  season.  Output  diagrams  are  written 
as  DISSPLA  metafiles. 

PC  Plotting.  Output  from  the  ocean  simulation 
model  consists  of  binary  and  ASCII  files.  To  avoid 
problems  associated  with  transferring  files  between 
different  hardware  architectures,  binary  files  should 
be  converted  into  ASCII  format  before  being  trans¬ 
ferred  to  PCs.  However,  ASCII-formatted  files  require 
more  disk  storage  than  binary.  If  storage  constraints 
prohibit  binary-to- ASCII  conversion,  commercial  and 
public  domain  software  is  available  (e.g.,  Kermit)  for 
binary  file  transfers  between  different  computers.  The 
export  programs  described  earlier  facilitate  data 
transfer  to  PCs. 

B.  Program  Structure 

Ehagrams  that  show  program  structure  and  provide 
comprehensive  details  of  the  main  program  and  its 
associated  subprograms  are  presented  in  the  appendix. 
Also  included  are  program  layering  diagrams,  which 
provide  a  summary  of  the  level  of  external  subroutine 
calls.  These  diagrams  provide  the  programmer  with 
an  outline  of  the  manner  in  which  the  subprograms 
are  called,  the  input  and  output  functions,  and  the 
logic  flow. 

C.  Data  Structures 

Data  structures  exist  as  binary,  ASCII,  and 
DISSPLA  metaplot  files.  These  files  are  used  as 
program  input,  are  created  as  program  output,  and 
are  stored  as  disk  files  on  the  resident  computer 
system.  Interactive  program  prompts  acquire  necessary 
information,  other  than  that  stored  as  disk  files,  for 
successful  program  execution. 

Input  Files.  LEVITUS.DAT,  a  climatological  atlas 
of  the  world  ocean  is  based  on  Levitus  (1982).  The 


ocean  simulation  model  requires  LEVITUS.DAT  as 
input.  LEVITUS.DAT  is  an  unformatted,  direct 
access  file  requiring  7.29  Mbytes,  or  14,580  VMS 
blocks,  of  disk  storage.  The  record  structure  consists 
of  180  fields  of  REAL*4  data,  including.numbcr  of 
temperature  observations,  mean  temperature  (in  units 
of  degrees  centigrade),  standard  deviation,  number  of 
salinity  observations,  mean  salinity  (in  units  of  parts 
per  thousand  or  practical  salinity  units),  and  standard 
deviation. 

Output  Files.  MODELl.DAT  is  an  unformatted 
direct  access  file  that  contains  displacement,  modified 
temperature  and  salinity,  and  sound  velocity  fields. 
The  record  structure  consists  of  4*NZ  fields  of 
floating  point  data  where  NZ  is  the  number  of  depths. 
Table  2  summarizes  the  output  files  for  this  model. 

MODEL l.AUX  is  an  unformatted,  direct  access 
file  that  contains  defining  parameters  fiabeled)  that 
should  be  used  to  augment  data  in  MODELl.DAT. 
The  record  structure  consists  of  3*NZ  fields,  where 
NZ  is  the  number  of  depths. 

MODEL l.UV  is  an  unformatted,  direct  access  file 
that  contains  the  velocity  fields  U,  V,  and  W.  TTie 


record  structure  contains  3*NZ  fields,  where  NZ  is 
the  number  of  depths. 

MODEL l.EIG  is  an  unformatted,  direa  access  file 
that  contains  modal  eigenvalues  and  eigenfunctions 
for  W  and  k.  The  record  structure  contains  NZ+1 
fields,  where  NZ  is  the  number  of  depths. 

MODELl.CTL  is  an  unformatted,  direct  access  file 
that  contains  control  information  relating  to  program 
restart  and  eigenmode/eigenvalue  control  and  use.  The 
record  structure  contains  five  fields. 

DISSPLA  META-PLOT  FILES  contain  a  non- 
device-specific  file  of  data  from  which  plots  can  be 
generated  using  a  device-specific  postprocessor.  The 
software  used  for  terminal  and  hardcopy  graphics  is 
Computer  Associates’  DISSPLA.  The  intermediate 
meiaplot  files  are  named  POPFIL.DAT. 

V.  Usage 

The  ocean  simulation  model  has  been  executed  on 
DEC  VAX/VMS  11/750  and  MicroVAX  11  systems 
and  on  the  CONVEX  UNIX  system.  Adherence  to 
most  Fortran  77  standards  has  enabled  program 
portability  between  the  VAX  and  CONVEX  systems. 


Table  2.  Ocean  Simulation  Model  Output  tiles. 


file  name 

FMT 

ACC 

RECL 

VARIABLES 

RECORD  STRUCTURE 

MODELl.DAT 

U 

D 

4*NZ 

Real 

U.V.W. 

Z 

(U(i),V(i),W(i), 

Z(i),  i.I.NZ) 

MODEL1.AUX 

U 

S 

3*NZ 

Real 

NT.NX.NZ.DT. 

DX,DZ,T0,LAT, 

LON.AZ 

NT.NX.NZ.DT, 

DX.DZ.TO.LAT, 

LON.AZ 

MODEL1.UV 

U 

D 

3’NZ 

Real 

u.v.w 

i-1,NZ) 

(U(i),V(i),W{i). 

MODEL1.EIG 

U 

D 

NZ+1 

Real 

k  W 

k.(W(i),i=1,NZ) 

MODELl.CTL 

U 

D 

5 

Rec  1, 

Int  NX 

Real  DX.XMAX 

NX, DX.XMAX 

Rec  2, 

Int  NZ 

Real  DZ.ZMAX 

NZ.DZ.ZMAX 

Rec  3, 

Int  NT 

Real  DT.TMAX 

NT.DT.TMAX 

Rec  4, 

Int  LAT.LON 

LATI.LONI 

Real  AZIMUTH 

LAT.LON. LATI.LONI, 

AZIMUTH 

Rec  5, 

Int.  ITO.NDIR 

Real  T 

T.ITO.NDIR 

Notes: 

-  File  format  (FMT)  is  denoted  as  unformatted.  U 

-  File  access  (ACC)  is  denoted  as  D  for  direct  access,  or  S  for  Sequential 
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A.  VMS  (DEC-VAX) 

The  ocean  simulation  model  can  be  executed  either 
interactively  or  in  batch  mode.  The  Levitus  climato¬ 
logical  atlas,  LEV1TUS.DAT,  is  required  as  part  of 
program  input.  Interactive  prompts  acquire  remaining 
necessary  program  information. 

Execution  of  the  MODELl  simulation  program 
with  the  input  described  in  section  V.E  on  the 
VAX  11/750  required  47.49  minutes  of  processing 
time.  The  1 1/750  possessed  7  Mbytes  of  memory  and 
a  floating-point  accelerator.  The  system  (SYSGEN) 
parameter.  VIRTUALPAGECOUNT,  was  set  to 
21,384.  User  accounts  required  modification  of 
the  parameter  PGFLQUO  from  16,000  blocks  to 
20,000  blocks. 

Because  the  11/750  was  used  by  other  interactive 
users,  executions  were  performed  in  batch  mode.  This 
method  made  it  possible  to  avoid  resource  depletion 
and  slowing  down  of  the  system  for  other  interactive 
users. 

On  the  Micro  VAX  II,  the  same  program  and  input 
required  32.60  minutes.  This  MicroVAX  is  a 
standard  system  containing  10  Mbytes  of  memory  and 
a  floating-point  accelerator.  The  SYSGEN  parameter, 
VIRTUALPAGECOUNT,  was  set  to  20,000.  Users’ 
accounts  required  setting  the  parameter  PGFLQUO 
from  16,000  blocks  to  20,000  blocks. 

B.  UNIX 

The  ocean  simulation  model  has  been  recompiled 
and  executed  on  a  CONVEX  computer  system.  The 
source  code  required  minimal  revision  to  reconcile 
the  record  length  parameters,  as  discussed.  The  source 
code  was  compiled  using  the  following  command  line 

fc  -02  model l.f  -o  model  1. 

The  test  program  and  input  ran  in  3.38  minutes. 
To  run  the  model,  the  associated  files  had  to  be  deleted 
before  starting  the  model.  For  details,  see  the  examples 
in  section  V.  E. 

C.  MS-DOS 

The  ocean  simulation  model  has  not  been  tested 
on  a  PC.  The  degree  of  software  modification  is 
dependent  upon  the  PC  Fortran  implementation. 
Model  output  can  be  converted  to  ASCII  format  and 
down  loaded  onto  PCs  for  plotting  and  listing. 

D.  Program  Input  Prompt  Descriptions 

The  ocean  simulation  model  requires  user  input 
for  execution.  If  MODELl  has  already  been  executed. 


then  existing  output  files  are  used  to  avoid  time- 
consuming  recalculations.  The  file  MODEL l.EIG 
contains  eigenvalue/eigenfunction  data  and  can  be 
reused.  Restart  information  is  saved  in  MODELl. CTL. 
MODELl  automatically  searches  for  the  existence  of 
these  files  before  computation.  The  following  input 
is  required. 

Latitude  and  Longitude  Prompts:  Floating-point 
values  for  latitude,  longitude,  and  Erection  of  section 
arc  required.  Invalid  coordinates  are  flagged  and  the 
prompt  is  reissued. 

Season  Prompt:  A  character  string  is  required.  The 
valid  responses  are  WINTER,  SPRING,  SUMMER 
or  FALL.  Because  MODELl  prompt  responses  are 
case  sensitive,  only  upper-case  values  are  vahd. 

Maximum  and  Delta  Ranges  in  km  Prompt:  The 
range  is  the  horizontal  distance  (in  kilometers)  of 
the  section.  The  maximum  range  is  the  total  length 
of  the  section,  and  delta  range  provides  the  horizon¬ 
tal  resolution  for  the  calculation.  In  the  present  version 
of  MODELl,  the  total  number  of  ranges  must  be  less 
than  or  equal  to  512. 

Maximum  and  Delta  Depths  in  m  Prompt:  The 
maximum  and  delta  depths  correspond  to  ranges  but 
apply  to  the  vertical  coordinate.  If  the  depth  incre¬ 
ment  is  chosen  such  that  the  number  of  depth 
intervals  would  be  greater  than  512,  then  the  pro¬ 
gram  sets  the  number  of  intervals  to  512  and  rescales 
the  depth  increment  accordingly. 

Maximum  and  Delta  Times  in  s  Prompt:  The 
program  allows  the  time  series  of  the  vertical  sec¬ 
tions  to  be  calculated.  The  times  are  specified  in 
seconds.  No  limits  are  placed  on  the  ranges  of  the 
input  times,  but  if  the  number  of  time  increments  is 
too  large,  then  the  available  storage  on  the  computer 
system  can  be  exceeded. 

Number  of  Directions  in  Isotropic  Case  Prompt: 
For  the  isotropic  case  (the  only  one  presently  allowed), 
each  mode  and  frequency  can  have  up  to  50  plane 
waves  propagating  at  different,  randomly  chosen 
directions.  The  number  of  such  waves  is  a  user  input. 
For  most  purposes,  NDIR  =  3  seems  to  produce 
reasonable  results. 

Number  of  Eigenvalues  (NEIG),  Modes 
(NMODES),  and  Frequencies  (NF)  prompt:  Inte¬ 
ger  values  for  NEIG,  NMODES,  and  NF  arc  required. 
Default  values  will  be  used  if  zeroes  are  entered  where 
NEIG  =  1000,  NMODES  =  5,  and  NF  =  8.  If  NEIG 
is  not  set  to  I,  MODELl  will  compute  as  many 
eigenfunctions  per  mode  and  frequency  as  there  are 
steps  in  the  horizontal.  If  NEIG  =  1,  only  the  first 
position  will  be  used  to  determine  the  Brunt- VaisiUa 
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frequency  profile  and  associated  eigenfunctions.  If 
NMODES  =  0  and  NF  =  0,  the  defaults  will  be 
substituted;  otherwise,  NMODES  and  NF  will  equal 
the  user-specified  value. 

Garrett-Munk  Brunt-V^s^a  Frequency  Profile 
Prompt:  A  character  string  (YES  or  NO)  is  required. 
If  YES,  a  generic  Brunt-VaisSia  frequency  profile  is 
used  in  place  of  the  real  profile. 

M2  Internal  Tide  Prompt:  A  character  string 
(YES  or  NO)  is  required.  If  YES,  a  generic  internal 
tide  (M2)  is  included  in  the  model.  In  this  case,  a 
1-m  internal  tidal  amplitude  will  be  assumed. 

At  the  end  of  input  prompts,  MODEL  1  displays 
some  of  the  user-supplied  information  and  computed 
values. 

E.  Sample  Program  Execution 
VAXA^MS.  a  sample  program  execution  with 
inputs  follows  where  program  input  is  shown  in 
boldface  type.  Additional  run-time  information  is 
available  in  the  General  Description  Section  for 
MODEL  1.  The  program  produces  some  output  other 
than  the  prompts.  This  output  is  intended  to  guide  the 
user.  For  example,  if  the  RESTART  parameter  is  false, 
then  it  tells  the  user  that  the  program  will  be  reini¬ 
tialized.  If  it  is  true,  then  the  user  has  the  choice  of 
reinitializing  the  program  or  continuing  from  where 
it  left  off.  The  •EXIST  parameter  tells  the  user 
whether  or  not  the  file  continues  to  exist.  If  only 
some  exist,  it  is  probably  be  better  to  delete  all  the 
files  and  reinitialize.  After  all  inputs  are  satisfied,  a 
summary  of  the  parameters  that  will  be  used  to 
generate  the  model  are  listed.  These  parameters  are 
for  information  only  and  play  no  part  in  the  subsequent 
operation  of  the  model. 

$  RUN  MODELl 

RESTART  =  F 
AUXEXIST  F 
DATEXIST  F 
CTLEXIST  F 
EIGEXIST  F 
RESTART  =  F 

••••♦••OCEAN  SIMULATION  MODEL^^^^^^^ 
VERSION  1.0 

Enter  latitude,  longitude  and  direction  of  section  in 
decimal  degrees. 

30,-70,123 

Enter  season  (WINTER,  SPRING,  SUMMER  OR 
FALL) 


SUMMER 

Enter  maximum  and  delta  ranges  in  km, 
maximum  and  delta  depths  in  m, 
maximum  and  delta  times  in  s 

40,.l 

3000,10 

100^00. 

ENTER  THE  NUMBER  OF  DIRECTIONS  TO  USE 
IN  ISOTROPIC  CASE 
3 

ENTER  NEIG,  NMODES,  NF 

1,0,0 

DO  YOU  WANT  GENERIC  GARRET-MUNK  BVF 
PROFILE? 

NO 

DO  YOU  WANT  M2  INTERNAL  TIDES? 

NO 

INPUT  DATA 

Latitude  =  30.000 
Longitude  =  -70.000 
Azimuth  =  123.000 

Xmax  =  40.000 

Zmax  =  3000.000 

Tmax  =  100. 

NX,NZ,NTJ^DIR  =  401  301  1  1 

COMPUTED  VALUES 

LAT  =30.000  LON  =-70.000 
LATl  =  25.000  LONl  =  -65.000 
Reduced  Azimuth  =  135.000 

Output  from  this  simulation  can  be  reviewed  with 
the  diagnostic  programs.  Figure  1  shows  a  VMS 
command  file  that  executes  the  program  MODELl. 

CONVEX/UNIX.  Example  shell  files  listings  to 
run  the  model  are  shown  in  Figure  2.  These  files  are 
mutually  interdependent.  One  important  note  in 
running  the  program  on  the  CONVEXAJNIX  system 
is  that  the  LEVITUS.DAT  file  must  be  named  in 
upper-case  characters  and  must  reside  in  the  same 
directory  (subdirectory)  as  the  MODELl  program. 

F.  General  Description 
MODELl.  The  ocean  simulation  model  is  an 
interactive,  computation-intensive  program.  Because 
of  the  high  number  of  iterations  needed  by  eigen¬ 
function/eigenvalue  computations,  MODELl  is 
computer  bound.  To  save  processing  time,  MODELl 
can  use  existing  files  from  previous  executions. 
Significant  time  can  be  saved  by  reusing  output  files. 
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$!  VAX/VMS  COMMAND  FILE  FOR  EXECUTING 
$!  OCEAN  SIMULATION  MODEL,  MODEL1 
$! 

$!  REMOVE  ANY  PRE-EXISTING  OUTPUT  FILES 
SPURGE 

SDEL  DEBUG.DAT;* 

SDEL  DIAGN0ST1CS.DAT;* 

SDEL  MODEL1.DAT:* 

SDEL  MODEL1.AUX;* 

SDEL  MODEL1.EIG;* 

SDEL  MODEL1.CTL;' 

SDEL  MODEL1.UV:* 

S! 

S! 

S!  BEGIN  EXECUTION 
SRUN  MODEL1 
30,-70,135 
SUMMER 
40,.  1 
3000,10 
1,100 
3 

1,0,0 

NO 

NO 


Figure  1.  VAX/VMS  Command  file  listing  for  running 
the  test  cose  for  MODElI  . 


File:  modrun 
# 

clean 

modell  <  input.dat 
time 

File:  clean 

#  cd  \bigone\briggs 
Is 

rm  MODELI. CTL 
rm  MODEL1.DAT 
rm  MODELI. AUX 
rm  MODELI. EIG 
rm  MODELI  . UV 
rm  DEBUG.DAT 
rm  DIAGNOSTICS.LIS 
Is 

File:  input.dat 

30,-70,135. 
SUMMER 
40, 10 
3000,10. 

1,100. 

3 

1,0,0 

NO 

NO 


Figure  2.  CONVEXiUNIX  Shell  command  file  listings 
for  running  the  test  case  for  MODELI. 


Automatic  inquiry  is  made  for  the  presence  of 
MODELI  output  files.  If  all  files  are  found  to  exist, 
then  MODELI  will  attempt  to  use  them.  An  error 
will  occur  if  these  files  were  produced  from  an  execution 
with  an  unsuccessful  program  terminatirai.  When  all  files 
do  not  previously  exist,  MODELI  creates  a  new  set. 

Section  V.E  contains  a  sample  program  execution. 
Note  that  immediately  after  starting  the  program, 
MODELI  reports  whether  previous  output  files  exist. 
An  F  (FALSE)  indicates  no  existence  and  T  (TRUE) 
indicates  the  files  exist. 

Storage  requirements  increase  as  values  for  the 
number  of  modes  (NMODES)  and  the  number  of 
frequencies  (NF)  increase.  Array  elements  are  computed 
with  the  formula  NMODES*NF*NZ*NX,  where  NZ  is 
the  number  of  depths  and  NX  is  the  number  of  points 
on  the  horizontal  region. 

MODELl_EXPORT.  Initial  information  about 
time  (NT),  horizontal  region  (NX),  depth  (NZ),  and  delta 
values  (DT,  DX,  DZ),  are  displayed.  Prompts  are  used 
to  acquire  desired  plot  type  and  iteration  number.  These 
inputs  are  used  as  indices  for  reading  MODELI  .DAT. 
Output  is  written  in  ASCII  format  to  a  file  named 
MODELI. PC. 

MODELl_SV_EXPORT.  Only  one  prompt,  for  iter¬ 
ation  number,  is  issued  Data  ccxitained  in  MODELl.DAT 
are  used  to  compute  sound  velocity.  Output  consisting 
of  horizontal,  depth,  and  sound  velocity  data  is  written 
to  an  ASCQ-formatted  file  named  MODELl_SV_ASCn. 

LEVASC.  No  interactive  input  is  required.  LEVASC 
reads  LEVITUS.DAT  (direct  access)  and  writes  out  the 
contents  in  ASCII  format  to  a  file  named  LEVTTUS.ASC. 

ASCLEV.  No  interactive  input  is  required.  ASCLEV 
reads  LEVTTUS.ASC  (ASCII  format)  and  writes  out 
the  contents  in  direct  access  format  to  a  file  named 
LEVITUS.DAT. 

MODELl_PLOT.  High-resolution  color  plots  are 
available  on  the  DEC  GKS  workstation. 
MODELl_PLOT  is  an  interactive  program  that 
prompts  for  desired  fields  and,  using  the  file 
MODELl.DAT,  computes  plot  fields  for  display  on  the 
GKS  high-resolution  color  monitor. 

MODELl_CONTOUR.  Prompts  are  used  to  acquire 
the  name  of  the  file  and  fields  to  be  plotted.  DISSPLA 
metafiles  that  contain  contour  plots  are  generated  by 
MODELl_CONTOUR.  Figure  3  contains  several 
examples  of  program  execution.  The  output  metafile 
must  be  postprocessed  for  device  specifications.  Sample 
output  plots  are  displayed  in  Figures  4  through  7. 

MODELl_ENERGY.  Only  one  prompt  is  issued. 
MODELl_ENERGY  automatically  opens  files 
(MODELI. AUX,  MODELl.DAT  and  MODELI. UV) 
that  contain  data  used  in  computations.  If  a  generic 
Garrett-Munk  Brunt-VUisHU  profile  is  used,  then  a  reply 
of  YES  should  be  used.  (In  all  other  cases,  a  reply  of 
NO  should  be  used.)  In  this  case,  the  computed  energy 
profiles  are  estimated.  Summary  information  is  written 
to  an  ASCn  file  named  ENERGY. US. 


$!  COMMAND  FILE  FOR  USING  MODEL1_CONTOUR  TO  CREATE  PLOT 
$!  FILES  OF  OCEAN  SIMULATION  MODEL1  FIELDS  ZO.  TO,  SO,  AND  SV 
$! 

$!  NOTES; 

$!  “STD00001.DAT"  IS  THE  DISSPLA  DEVICE  SPECIFIC  FILE 
$!  NAME  FOR  THE  DEC  LN03  LASER  PRINTER. 

$' 

$!  DISSPLA  APPENDS  PLOTS  TO  EXISTING  DEVICE  SPECIFIC  FILES, 

$!  DELETE  OR  RENAME  EXISTING  VERSIONS  TO  AVOID  UNINTENTIONAL 
$!  PLOT  APPENDING. 

$! 

$!  "DISSIOS"  IS  A  COMMAND  WHICH  INVOKES  THE  DISSPLA 
$1  POST-PROCESSOR.  THE  POST  PROCESSOR  GENERATES 
$!  DEVICE  SPECIFIC  PLOT  FILES.  HERE,  THE  LN03  SPECIFIC 
$!  PLOT  FILE  IS  NAMED  STD00001.DAT 
$! 

$!  "LNOS"  IS  A  COMMAND  THE  SENDS  A  FILE  TO  THE  LN03 
$!  PRINT  OUEUE 
$' 

$SET  VERIFY 

SDELETE  STD00001.DAT;- 
SPURGE 
$! 

SRUN  MODEL1_CONTOUR 
MODEL1 
ZD 

SDISSIOS 
12 

DRAW=1-END$ 

$LN03  STD00001.DAT 
S' 

SRUN  MODEL1_CONTOUR 
MODEL1 
TD 

SDELETE  STD00001.DAT:* 

SDISSIOS 
12 

DRAW=1-ENDSSRUN  MODEL1_CONTOUR 
SLN03  STD00001.0AT 
S' 

SRUN  MODEL1_CONTOUR  !  PLOT  SALINITY  ANOMALY 

MODEL1 

SD 

SDELETE  STD00001.DAT:* 

SDISSIOS 

12 

DRAW =1 -ENDS 
SLN03  STD00001.DAT 
SI 

SRUN  MODEL1_CONTOUR  '  PLOT  SOUND  VELOCITY  ANOMALY 

MODEL1 

SV 

SDELETE  STD00001.DAT;* 

SDISSIOS 

12 

DRAW  =1 -ENDS 
SLN03  STD00001.DAT 
SPURGE 


Figure  3.  VAX/VMS  Command  file  for  producing  contour  plots  of fields  generated  by  MODELl . 


!  AVOID  PLOT  APPENDING 


!  PLOT  VERTICAL  DISPLACEMENT 
!  INPUT  DATA  FILE  NAME 
I  FIELD  TO  PLOT 
!  CREATE  DEVICE  SPECIFIC  FILE 
!  POST-PROCESSOR  DIRECTIVE 
!  POST-PROCESSOR  DIRECTIVE 
!  PRINT  THE  PLOT  FILE 

!  PLOT  TEMPERATURE  ANOMALY 
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Figure  4  Contour  plot  of 


Figure  5  Contour  plot  of  h 


X-DISTANCE  (km) 


Figure  6.  Contour  plot  of  salinity  anomaly,  produced  by  MODELl  CONTOUR.  The  ragged  contour  at  the  bottom  of  the  plot  is  an 
artifact  produced  by  a  nearly  zero  anomaly.  Contour  interval  =  0.005. 


X-DISTANCE  (km) 


Figure  7.  Curuuur  plot  oj  sound  veiocity  anomaly,  produced  by  MODEL/  CONTOUR.  A  slowly  varying  field  near  zero  causes  the 
ragged  contours  in  the  lower  part  of  the  diagram.  Contour  interval  =  0.200. 
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DISSPLA  metafiles  that  contain  internal  wave 
distribution  plots  are  generated.  Figure  8  shows  a 
sample  VAX/VMS  command  file  for  program 
execution.  Figure  9  shows  a  sample  plot  of  horizontal 
distribution.  Figure  10  shows  a  sample  plot  of 
vertical  distribution. 

MODELl_UVW_READ.  Time  (NT),  horizontal 
region  (NX),  depth  (NZ),  and  delta  values  (DT,  DX.DZ), 
are  displayed  first.  Prompts  are  used  to  acquire  infor¬ 
mation  about  decimation  rates  for  the  horizontal  and 
vertical  direction.  These  inputs  are  used  as  indices  for 
reading  U.V,  and  W  fields  in  MODEL  1. DAT.  Output  is 
written  to  the  user’s  terminal. 

MODELl_CTL_LIST.  There  are  no  input  prompts 
for  MODELl_CTL_LIST.  The  contents  of  the  control  data 
file,  MODELl.CTL,  are  displayed  on  the  user’s  terminal. 

MODELl_LOOK.  Input  prompts  are  used  to  acquire 
information  about  field  minimum,  maximum,  and  delta 
values.  These  values  are  used  as  indices  for  depth  and 
for  the  horizontal  region.  Model  displacement  fields  that 
contain  the  file  MODELl.DAT  arc  displayed  on  the 
user’s  terminal. 


MODELl_EIG_LOOK.  Input  prompts  acquire 
information  about  desired  iteration  number,  horizontal 
field,  number  of  frequencies,  and  number  of  modes. 
Plots  of  the  horizontal  field  are  generated  using  the 
desired  iteration  number  and  designated  frequencies  and 
modes.  Values  for  horizontal  fields  and  associated 
frequency,  mode,  and  wave  number  arc  listed. 

lVIODESUB_TEST.  Input  prompts  are  used  to 
acquire  information  about  position,  season,  depth, 
and  frequency.  Internal  wave  modes  for  Brunt-Vaisala 
frequency  profiles  are  computed  and  written  as  a 
DISSPLA  metafile.  Figure  1 1  contains  a  sample 
VAX/VMS  command  file  for  program  execution. 
Figure  12  contains  a  sample  output  plot. 

DISPERSION.  Prompts  are  used  to  acquire 
information  on  position,  season,  depth,  number  of 
frequencies,  and  number  of  modes.  Modal  dispersion 
diagrams  for  a  Brunt-Vaisala  frequency  profile  at  a  given 
location  and  season  arc  written  as  a  DISSPLA  metafile. 
Figure  13  contains  a  sample  VAX/VMS  command  file 
for  program  execution.  Figure  14  contains  a  sample 
output  plot. 


$'  COMMAND  FILE  FOR  USING  M0DEL1_ENER6Y  TO 


$ 

S' 

$f 
$! 

$' 

$' 

$' 

S' 

s 
s 
s 
s 
s 
s 

S' 

S' 

S' 

s 
s 
s 

S' 

SDELETE  STD00001.DAT;* 

S! 

SRUN  MODEL  1._ENERGY 

YES 

S' 

S' 

SDISSIOS 

12 

MODI  =  1  -END(HWR0T  =  M0VIE)S 
LN03  STD00001  DAT 


CREATE  PLOT  FILES  OF  OCEAN  SIMULATION  M0DEL1 
INTERNAL  WAVE  ENERGIES 

NOTES, 

•'ST000001.DAT  '  IS  THE  DISSPLA  DEVICE  SPECIFIC  FILE 
NAME  FOR  THE  DEC  LN03  LASER  PRINTER 

DISSPLA  APPENDS  PLOTS  TO  EXISTING  DEVICE  SPECIFIC  FILES 
DELETE  OR  RENAME  EXISTING  VERSIONS  TO  AVOID  UNINTENTIONAL 
PLOT  APPENDING 

"DISS105  '  IS  A  COMMAND  WHICH  INVOKES  THE  DISSPLA 
POST-PROCESSOR  THE  POST  PROCESSOR  GENERATES 
DEVICE  SPECIFIC  PLOT  FILES.  HERE.  THE  LN03  SPECIFIC 
PLOT  FILE  IS  NAMED  STD00001.DAT 

"LN03'-  IS  A  COMMAND  THE  SENDS  A  FILE  TO  THE  LN03 
PRINT  QUEUE 


'  AVOID  PLOT  APPENDING 


'  PROGRAM  INPUT 


A  GENERIC  BRUNT- 
VAISALA  PROFILE 


'  CREATE  DEVICE  SPECIFIC  FILE 
'  POST-PROCESSOR  DIRECTIVE 
'  POST  PROCESSOR  DIRECTIVE 
'  PRINT  THE  PLOT  FILES 


Figure  8  VAXiVMS  Command  file  listing  for  running  MODELl  ENERGY  and  plotting 
the  results. 


$!  COMMAND  FILE  FOR  USING  MODESUB_TEST  TO 
$!  CREATE  A  PLOT  FILE  OF  OCEAN  SIMULATION  MODEL1 
$!  INTERNAL  WAVE  MODES  FOR  BRUNT-VAISALA  FREQUENCY 
$!  PROFILES 
$! 

$!  NOTES: 

$!  "STD00001.DAT"  IS  THE  DISSPLA  DEVICE  SPECIFIC  FILE 
$!  NAME  FOR  THE  DEC  LN03  LASER  PRINTER, 

$! 

$!  DISSPLA  APPENDS  PLOTS  TO  EXISTING  DEVICE  SPECIFIC  FILES. 

$!  DELETE  OR  RENAME  EXISTING  VERSIONS  TO  AVOID  UNINTENTIONAL 
$!  PLOT  APPENDING. 

$! 

$!  "DISSIOS”  IS  A  COMMAND  WHICH  INVOKES  THE  DISSPLA 
$!  POST-PROCESSOR.  THE  POST  PROCESSOR  GENERATES 
$1  DEVICE  SPECIFIC  PLOT  FILES.  HERE,  THE  LN03  SPECIFIC 
$1  PLOT  FILE  IS  NAMED  STD00001.DAT 
$! 

$!  INOS"  IS  A  COMMAND  THE  SENDS  A  FILE  TO  THE  LN03 
$!  PRINT  QUEUE 
$! 

$! 

$! 

SDELETE  STD00001.DAT;* 

$' 

$R  MODESUB_TEST 
30.-70 
SUMMER 
3000.10 
2.0 
$' 

SDISSIOS 
12 

M0DI  =  1-END(HWR0T=M0VIE)$  !  POST-PROCESSOR  DIRECTIVE 
LN03  STD00001  OAT  !  PRINT  THE  PLOT  FILE 


!  AVOID  PLOT  APPENDING 


'  PROGRAM  INPUT  -  LOCATION 
!  PROGRAM  INPUT  -  SEASON 
!  PROGRAM  INPUT  -  DEPTH  MAX  &  MIN 
!  PROGRAM  INPUT  -  FREQUENCY 

!  CREATE  DEVICE  SPECIFIC  FILE 
!  POST-PROCESSOR  DIRECTIVE 


Figure  II.  VAX/VMS  Command  file  for  running  MODESUB  _TEST. 
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-1.5  0  1.5 

MODE  AMPLITUDE 
FREQUENCY  (cph)  2.000 


0  5  10  15 

FREQUENCY 

Figure  12.  Plot  output  from  the  command  file  in  Figure  11. 
showing  the  modes  for  the  Bnmt-Vdisdld  frequency  profile 
at  30°N,  70°W  at  2  cph,  which  is  plotted  on  the  right. 
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COMMAND  FILE  FOR  USING  DISPERSION  TO 
CREATE  A  PLOT  FILE  OF  OCEAN  SIMULATION  MODEL1 
MODAL  DISPERSION  DIAGRAMS 

NOTES: 

''STD00001.DAT”  IS  THE  DISSPLA  DEVICE  SPECIFIC  FILE 
NAME  FOR  THE  DEC  LN03  LASER  PRINTER. 

DISSPLA  APPENDS  PLOTS  TO  EXISTING  DEVICE  SPECIFIC  FILES. 
DELETE  OR  RENAME  EXISTING  VERSIONS  TO  AVOID  UNINTENTIONAL 
PLOT  APPENDING. 

■‘DISS105''  IS  A  COMMAND  WHICH  INVOKES  THE  DISSPLA 
POST-PROCESSOR.  THE  POST  PROCESSOR  GENERATES 
DEVICE  SPECIFIC  PLOT  FILES.  HERE,  THE  LN03  SPECIFIC 
PLOT  FILE  IS  NAMED  STDOOOGI.DAT 

“LN03”  IS  A  COMMAND  THE  SENDS  A  FILE  TO  THE  LN03 
PRINT  QUEUE 


SDELETE  STD00001.DAT;* 

$! 

$RUN  DISPERSION 
30,-70 
SUMMER 
3000,10 
100 
5 
$! 

$DISS105 

12 

MODI = 1  -ENO(HWROT= MOVIE)$ 
LN03  STD00001.DAT 


!  AVOID  PLOT  APPENDING 


!  PROGRAM  INPUT  -  LOCATION 
!  PROGRAM  INPUT  -  SEASON 
!  PROGRAM  INPUT  -  DEPTH  MAX  &  MIN 
!  PROGRAM  INPUT  -  #  OF  FREQUENCIES 
!  PROGRAM  INPUT  -  #  OF  MODES 

!  CREATE  DEVICE  SPECIFIC  FILE 
!  POST-PROCESSOR  DIRECTIVE 
!  POST-PROCESSOR  DIRECTIVE 
!  PRINT  THE  PLOT  FILE 


Figure  13.  VAXIVMS  Command  file  listing  for  running  DISPERSION. 
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0  12  3  4 

HORIZONTAL  WAVENUMBER  (cp  km) 


FREQUENCY 


Figur-  <4.  Plot  of  the  output  from  DISPERSION,  showing 
the  tine  tr  internal  wave  dispersion  diagram  for  the  Brunt- 
Vdisdld  frequency  profile  plotted  on  the  right. 


21 


n 


VI.  Conclusions 

The  first-level  ocean  simulation  was  designed  to 
provide  a  capability  for  simulating  the  perturbations 
to  the  temperature,  salinity,  and  velocity  fields  due  to 
internal  wave  and  tidal  activity  in  midocean  condi¬ 
tions.  It  was  designed  to  be  portable  and  easy  to  use 
and  maintain.  The  simulation  was  therefore  written  in 
standard  FORTRAN??  and  has  been  tested  on  Digital 
Equipment  Corporation’s  VAX  computers,  using  the 
VMS  operating  system,  and  on  a  CONVEX  computer, 
using  the  UNIX  operating  system.  The  program  should 
operate  properly  on  any  other  system,  provided  the 
record  lengths  of  the  direct  access  files  are  suitable. 

VII.  Recommendations 

This  program  should  be  used  for  acoustic  field 
simulation,  cruise  planning,  mooring  motion  simula¬ 
tion,  and  general  ocean  system  modeling.  Further 
development  of  the  internal  tidal  stmeture,  particularly 
the  propagation  of  the  internal  tide  in  shallow  coastal 
regions,  is  needed  to  make  the  simulation  more 
realistic  and  useful  on  continental  shelf  regions.  The 
inclusion  of  varying  bottom  depths  and  variable  output 
data  structures  would  make  the  program  more  useful 
to  the  acoustic  modeling  community. 
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MODELl  Program  Structure 


MAIN  PROGRAM 

I 

CONTROL_INPUT 

OPEN  OUTPUT  FILES 
EIG_CONTROL 
I 

PROFILE_INPUT 

I 

INT_WAVE_SIMULATION 
(TIME  LOOP) 

(X  -  LOOP) 

DISPLACE 
PROFILE_CALC 
(Z  -  LOOP) 

(OUTPUT  — >  MODELl.DAT) 
(END  ZLOOP) 

(END  X  -  LOOP) 

(END  TIME  -  LOOP) 
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MODELl  Program  Layering 


CONTROL  INPUT 


PROFILE  INPUT 


INT  WAVE  SIMULATION 


OPEN  OUTPUT  I 
FILES  I 

EIG  CONTROL 


INTRPL  BVFREQ  DIST 


DISPLACE 


PROFILE  CALC 


MODESUB 


INTERP 


II  II 

TURN  NUMEROV  INTERP  AVGINT 
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M0DEL1_EXP0RT  Program  Structure 


MAIN 

I 

PARSE  (MODELl.AUX) 

I 

INTERACTIVE  TERMINAL  INPUT 

I 

OPEN  MODELl.DAT  &  MODELl.UV 

I 

COMPUTE  PLOT  FIELDS 

I 

WRITE  FIELDS  =>  MODELl . PC 


MODELl_EXPORT  Program  Layering 


Main 


1  1  1 

PARSE  1  1 

1 

1 

1 

1 

MODELl.AUX  1  1 

1 

1 

INTERACTIVE  I 

1 

1 

INPUT  1 

1 

1 

OPEN  MODELl.DAT 

1 

1 

&  MODELl.UV 

1 

1 

COMPUTE  PLOT 

1 

FIELDS 

1 

WRITE  FIELDS 
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M0DEL1_SV_EXP0RT  Program  Structure 

MAIN 

I 

PARSE  (MODELl.AUX) 

I 

OPEN  MODELl.DAT 

I 

COMPUTE  SVEL  (SOUND  VELOCITY) 

I 

OUTPUT  =>  ASCII  FILE  CONTAINING 
SOUND  VELOCITY  FIELD 


MODELl_SV_EXPORT  Program  Layering 


MAIN 


PARSE 

1 

1 

MODELl.AUX 

1 

1 

OPEN 

COMPUTE 

MODELl.DAT 

SOUND  velocity 

OUTPUT  =>  ASCII  FILE 
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LEVASC  Program  Structure 


MAIN 

I 

OPEN  LEVITUS.DAT  (BINARY  FORMAT) 

I 

LOOP 

READ  LEVITUS.DAT 

I 

WRITE  LEVITUS.ASC  (ASCII  FORMAT) 
END  LOOP 


LEVASC  Program  Layering 

MAIN 

I 


OPEN  LEVITUS.DAT 

LOOP 

READ 

1 

LEVITUS.DAT 

1 

WRITE 

LEVITUS.ASC 
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ASCLEV  Program  Structure 


MAIN 

I 

OPEN  LEVITUS.ASC  (ASCII  FORMAT) 

I 

LOOP 

READ  LEVITUS.ASC 

I 

WRITE  LEVITUS.DAT  (BINARY  FORMAT) 
END  LOOP 


ASCLEV  Program  Layering 


MAIN 


OPEN  LEVITUS.ASC  LOOP 

READ  LEVITUS.ASC 


WRITE  LEVITUS.DAT 


M0DEL1_PL0T  Program  Structure 


MAIN  PROGRAM 

PARSE  (MODELl.AUX) 

I 

INTERACTIVE  INPUT 

I 

OPEN  MODEL.DAT  &  MODELl.UV 

I 

ESTABLISH  PLOTTING  FEATURES 
UI S  $CREATE_COLOR_MAP 
UIS$CREATE_DI SPLAY 
UI  S  $CREATE__WI  NDOW 
UIS$HSV_TO_RGB 
UIS$SET_COLOR 
UISSSET_WRITING_INDEX 
UIS$SET_LINE_WIDTH 
UISSPLOT 

I 

PERFORM  PLOT 

I 

ESTABLISH  COLOR  SCALE 
UIS$PLOT 
UIS$SET_FONT 
UIS$NEW_TEXT_LINE 
UIS$TEXT 

I 

CREATE  COLOR  SCALE 
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M0DEL1_PL0T  Program  Layering 


MAIN 


PARSE  III  I 

MODELl.AUX  III  I 

INTERACTIVE  |  I  I 

INPUT  1  I  I 

OPEN  I  I 

MODELl.DAT  |  I 

&  MODELl.UV  I  I 

I  I 

PLOT  FEATURES  I 

UIS$CREATE_COLOR_MAP  I 

UIS$CREATE_DI SPLAY  I 

UIS$CREATE_WINDO  I 

UIS$HSV_TO_RGB  I 

UIS$SET_COLOR  I 

UIS$SET_WRITING_INDEX  | 
UISSSET_LINE_WIDTH  I 

UISSPLOT  I 

I 

PERFORM  PLOT 


EST.  COLOR  SCALE 

UISSPLOT 

UIS$SET_FONT 

UI S  $NEW_TEXT_LINE 

UISSTEXT 


CREATE 
COLOR  SCALE 


M0DEL1_C0NT0UR  Program  Structure 


MAIN 

I 

INTERACTIVE  TERMINAL  INPUT 

I 

PARSE  (MODELl.AUX) 

1 

OPEN  MODELl.DAT  &  MODELl.UV 

I 

ESTABLISH  PLOT  PARAMETERS 

I 

COMPUTE  PLOT  FIELDS 
SVEL 

I 

GRAPHICS  SUBROUTINES 
COMPRS 
BCOMMON 
PAGE 
S CMP LX 
AREA2D 
HEAD IN 
YNAME 
XNAME 
GRAF 
FRAME 
CONMAK 
CONLIN 
CONTUR 
ENDPL 
DONEPL 
I 

OUTPUT  =>  DISSPLA  PLOT  METAFILE 
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M0DEL1_C0NT0UR  Program  Layering 


MAIN 


1  1  1 

INTERACTIVE  |  I 

1 

1 

1 

1 

1 

1 

INPUT  1  1 

1 

1 

1 

PARSE  1 

1 

1 

1 

MODELl.AUX  1 

1 

1 

1 

OPEN 

1 

1 

1 

MODELl.DAT 

1 

1 

1 

MODELl .UV 

1 

1 

1 

ESTABLISH  PLOT 

1 

1 

PARAMETERS 

1 

t 

COMPUTE 

PLOT  1 

FIELDS 

1 

1 

SVEL 

1 

GRAPHICS 

SUBROUTINES 

I 

COMPRS 

BCOMMON 

PAGE 

S CMP LX 

AREA2D 

HEAD IN 

YNAME 

XNAME 

GRAF 

FRAME 

CONMAK 

CONLIN 

CONTUR 

ENDPL 

DONEPL 


OUTPUT  =>  DISSPLA 
PLOT  METAFILE 


M0DEL1_ENERGY  Program  Structure 


MAIN 

I 

INTERACTIVE  TERMINAL  INPUT 

I 

PARSE  (MODELl.AUX) 

I 

OPEN  MODELl.DAT  &  MODELl . UV 

I 

ESTABLISH  PLOT  PARAMETERS 

I 

COMPUTE  PLOT  FIELDS 

I 

OUTPUT  =>  ENERGY. LIS 

I 

GRAPHICS  SUBROUTINES 
COMPRS 
1 

REPEAT  FOR  EACH  PLOT 
SETCLR 
PAGE 
AREA2D 
TRIPLX 
HEIGHT 
HEADIN 
YNAME 
XNAME 
GRAF 
LINES 
LEGLIN 
CONLIN 
SCLPIC 
CURVE 
DASH 
DOT 
LEGEND 
I 

OUTPUT  =>  DISSPLA  METAFILE 
RESET 
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M0DEL1_ENRGY  Program  Layering 


MAIN 


1  1  III 

INTERACTIVE  I  III 

1 

1 

1 

1 

TERMINAL  |  III 

1 

1 

INPUT  PARSE  1  1  1 

1 

1 

MODELl.AUX  1  1  1 

1 

1 

OPEN  1  1 

1 

1 

MODELl.DAT  |  | 

1 

1 

MODELl.UV  1  1 

1 

1 

EST.  PLOT  1 

1 

1 

PARAMETERS  | 

1 

1 

COMPUTE 

1 

1 

PLOT  FIELDS 

1 

OUTPUT=> 

ENERGY. LIS 

1 

1 

1 

GFIAPHICS 

COMPRS 

I 

REPEAT  FOR  E.ACH  PLOT 
SETCLR 
PAGE 
AREA2D 
TRIPLX 
HEIGHT 
HEADIN 
YNAME 
XNAME 
GRAF 
LINES 
LEGLIN 
CONLIN 
SCLPIC 
CURVE 
DASH 
DOT 

LEGEND 

ENDPL 

I 

OUTPUT=>REGIS  FILE 

I 

RESET 


DONEPL 


M0DEL1_UVW_READ  Program  Structure 


MAIN 

I 

PARSE  (MODELl.AUX) 

I 

INTERACTIVE  TERMINAL  INPUT 
I 

OPEN  MODELl.UV 
I 

READ  MODELl.UV 

I 

OUTPUT  =>  UV  FIELDS 


MODELl_UVW_READ  Program  Layering 


MAIN 


1 

PARSE 

1 

1 

1 

1 

1 

1 

1 

1 

MODELl .AUX 

1 

1 

1 

1 

INTERACTIVE 

1 

1 

1 

INPUT 

1 

OPEN 

MODELl.UV 

1 

1 

1 

READ 

MODELl . UV 

1 

1 

1 

1 

1 

OUTPUT 

UV  FIELDS 
TO  TERMINAL 
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M0DEL1_CTL_LIST  Program  Structure 


Main 

I 

OPEN  MODELl.CTL  FILE 

I 

OUTPUT  =>  CTL  INFORMATION 


MODELl_CTL_LIST  Program  Layering 


MAIN 


OPEN  I 

MODELl.CTL  | 

OUTPUT  =>  CTL 
INFORMATION  TO  TERMINAL 


MODELl_LOOK  Program  Structure 


MAIN 

I 

PARSE  (MODELl.AUX) 

I 

OPEN  M0DEL1.DAT 

I 

READ  (MODELl.DAT  FILE) 

I 

INTERACTIVE  TERMINAL  INPUT 

I 

OUTPUT  MODAL  DISPLACEMENTS 


MODELl_LOOK  Program  Layering 


MAIN 


PARSE  III  I 

MODELl.AUX  III  I 

OPEN  I  (  I 

MODELl.DAT  I  I  I 

READ  I  I 

MODELl.DAT  |  1 

INTERACTIVE  I 

TERMINAL  INPUT  I 

OUTPUT  =>  MODAL 
DISPLACEMENTS 
TO  TERMINAL 
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M0DEL1_EIGL00K  Program  Structure 


MAIN 

I 

PARSE  (MODELl.AUX) 

I 

OPEN  MODELl.EIG 

I 

GRAPHICS  SUBROUTINES 

COMPRS 

I 

REPEAT  FOR  EACH  PLOT 
SETCLR 
PAGE 
AREA2D 
FASHION 
HEIGHT 
HEADIN 
XNAME 
YNAME 
GRAF 
ME S SAG 
INTNO 
CURVE 
ENDPL 

I 

OUTPUT  =>  DISSPLA  METAFILE 


DONPL 


M0DEL1_EIGL00K  Program  Layering 

MAIN 


PARSE  I  I 

MODELl.AUX  I  I 

OPEN  I 

MODELl.EIG  I 

GRAPHICS  SUBROUTINES 

I 

COMPRS 

I 

REPEAT  FOR  EACH  PLOT 

I 

SETCLR 

PAGE 

AREA2D 

FASHION 

HEIGHT 

HEADIN 

XNAME 

YNAME 

GRAF 

MESSAG 

INTNO 

CURVE 

ENDPL 

I 

OUTPUT  =>  DISSPLA  METAFILE 

I 

DONPL 
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MODESUB_TEST  Program  Structure 


MAIN 

I 

CONTROL_INPUT 

OPEN  TESTJ51AGNOST1CS.1.IS 

I 

PROFILE  INPUT 
OPEN  LEVITUS.DAT 
LINT 
INTBRPL 
BVFREQ 

OUTPUT  ->  TEST  DIAGNOSTICS. LIS 

I 

INTJtAVE  SIKUIATION 
DISPLACE 
.  ODBSUB 

GRAPHICS  SUBPROGRAMS 
SET  PLOT  FEATURES 
(LOOP  FOR  EACH  INTERNAL  WAVE  MODE) 
OUTPUT  ->  DISSPLA  METAFILE 
(END  LOOP) 

BRUNT-VAISALA  frequence  profile  PLOTS 
OUTPUT  ->  DISSPLA  METAFILE 
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r,  0  1  fi. -i  0  J 


hUG  21  '3 


14:55 


F'GiiL  .  .JO 


MODESUB_TEST  Program  Layering 


CONTROL_INPUT  | 

1  I 

OPEN  TEST_DIAGNOSTICS . LIS  I 

PROFILE  INPUT 


OPEN  LINT 

LEVITUS.DAT 


I  I  I  I  I 

INTRPL  BVFREQ  DIST  OUTPUT  =>  I 

TEST  DIAGNOSTICS. LIS  I 


INT  WAVE  SIMULATION 


DISPLACE 


MODESUB 


GRAPHICS 

SUBPROGRAMS 


NUMEROV  INTERP 


AVGINT  I  I 

SET  PLOT  FEATURES  REPEAT  FOR  EACH 
I  INTERNAL  WAVE  MODE 

COMPRS  I 


FREQ 

PROFILE 


SETDEV 

PAGE 

PHYSOR 

AREA2D 

SERIF 

SHDCHR 

HEIGHT 

HEADIN 

XNAME 

YNAME 

GRAF 

FRAME 

GRACE 


CURVE  PHYSOR 

I  AREA2D 

lUTPUT  =>  SERIF 

OPFIL.DAT  SHDCHR 

HEIGHT 
HEADIN 
XNAME 
YNAME 
GRAF 
FRAME 
CURVE 
DASH 
CURVE 
ENDGR 
ENDPL 

OUTPUT  =>  POPFIL.DAT 
DONEPL 


OUTPUT  => 
POPFIL.DAT 
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DISPERSION  Program  Structure 


MAIN 

I 

CONTROL_INPUT 

OPEN  TEST_DIAGNOSTICS.LIS 

I 

PROFILE_INPUT 

OPEN  LEVITUS.DAT 

INTERPL 

BVFREQ 

OUTPUT  =>  TEST_DIAGNOSTICS. LIS 

I 

I NT_WAVE_S I MULAT I ON 
DISPLACE 
MODESUB 

GRAPHICS  SUBPROGRAMS 
SET  PLOT  FEATURES 

(LOOP  FOR  EACH  INTERNAL  WAVE  MODE) 
OUTPUT  =>  DISSPLA  METAFILE 
POPFIL.DAT 

(END  LOOP) 

BRUNT-VAISALA  FREQUENCY  PROFILE  PLOTS 
OUTPUT  =>  DISSPLA  METAFILE 


DISPERSION  Program  Layering 


MAIN 


I 

CONTROL  INPUT 


OPEN  TEST_DIAGNOSTICS.LIS  I 

PROFILE  INPUT 


INPUT  LINT  INTRPL  BVFREQ  DIST  OUTPUT  => 

LEVITUS.DAT  TEST  DIAGNOSTICS . LIS 


I 

I 


I NT  WAVE  SIMULATION 


DISPLACE 


MODE SUB 


TURN  NUMEROV  INTERP 


GRAPHICS 

SUBPROGRAMS 


AVGINT  I  I  I 

SET  PLOT  FEATURES  REPEAT  FOR  EACH  FREQ 

I  INTERNAL  WAVE  MODE  PROFILE 


COMPRS 

SETDEV 

PAGE 

PHYSOR 

AREA2D 

SERIF 

SHDCHR 

HEIGHT 

HEADIN 

XNAME 

YNAME 

GRAF 

FRAME 

GRACE 


CURVE  PHYSOR 

I  AREA2D 

OUTPUT  =>  SERIF 

POPFIL.DAT  SHDCHR 

HEIGHT 
HEADIN 
XNAME 
YNAME 
GRAF 
FRAME 
CURVE 
DASH 
CURVE 
ENDGR 
ENDPL 

OUTPUT  =>  POPFIL.DAT 
DONEPL 
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